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We consider a pair of three - level atoms interacting with a common vacuum and analyze the process of 
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collective damping can generate robust entanglement of the asymptotic states. 

SHORT TITLE: Stationary entanglement between three - level atoms 

PACS numbers: 03.67.Mn, 03.65.Yz, 42.50.-p 

Keywords: three - level atoms, collective damping, entanglement production 



I. INTRODUCTION 

The important problem of evolution of entanglement in re- 
alistic quantum systems interacting with their environments 
was mainly discussed in the case of two two - level systems 
(qubits). In that case, the interesting idea that dissipation can 
create rather then destroy entanglement was studied in details. 
In particular, in the case of two - level atoms, the possible 
production of robust or transient entanglement induced by the 
process of spontaneous emission was shown HI 0, S 0] • 

Much more complex and interesting is the process of cre- 
ation of entanglement involving multilevel atoms. In such a 
case, quantum interference between different radiative tran- 
sitions can influence the dynamics of the system. For a pair 
of largely separated three - level atoms the role of such in- 
terference in the process of degradation of entanglement was 
studied in Ref. |5]. When the interatomic distance is com- 
parable to the wavelength of the emitted radiation, the cou- 
pling between the atoms via common vacuum gives rise to 
the collective effects such as collective damping and dipole - 
dipole interaction. Such effects are well known |6], particu- 
larly in the case of two - level atoms. In the system of three - 
level atoms having closely lying excited states, radiative cou- 
pling can produce a new interference effect in the spontaneous 
emission. This effect manifests by the cross coupling between 
radiation transitions with orthogonal dipole moments [7]. All 
such collective properties of the system influence the entan- 
glement between three - level atoms. 

In the paper, we study the entanglement production be- 
tween three - level atoms due to the collective damping (a 
detailed analysis of the entanglement evolution in the pres- 
ence of all collective effects will be presented elsewhere). In 
that case, the analysis is involved since there is no simple nec- 
essary and sufficient condition of entanglement for a pair of d 
- level systems with d >3. Peres-Horodecki separability cri- 
terion [8, 9] only shows that states which are not positive after 
partial transposition (NPPT states) are entangled. But there 
can exist entangled states which are positive after this opera- 
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tion fl^ (bound entangled PPT states). The problem of exis- 
tence of bound entangled (i.e. non-distillable Lil]) states can 
be analyzed in terms of the rank of density matrix of bipartite 
system and the ranks of its partial traces. If a state is separable 
or bound entangled, then its rank must be larger then the ranks 
of paitial traces |12]. In the paper we focus on the possibility 
of creation of NPPT "free" entangled states, so we not discuss 
these problems. To detect and quantify entanglement we use 
negativity of partial transposition of density matrix. 

As we show, in the limit of small separation between the 
atoms, the process of the photon exchange between the atoms 
produces such correlations that the dynamics is not ergodic 
and there are nontrivial asymptotic stationary states. We com- 
pute the explicit form of the asymptotic state for any initial 
state and show that some of the asymptotic states are NPPT 
states, even if the initial states were PPT states. This effect 
occurs for example for a large class of diagonal i.e. separable 
initial states. We give also the example of bound entangled 
PPT state which evolves into NPPT (i.e. distillable) asymp- 
totic state. 



II. DYNAMICAL EVOLUTION OF TWO THREE-LEVEL 
ATOMS 

Consider two identical three - level atoms ( A and B) in 
the V configuration. The atoms have two near - degenerate 
excited states |2^) (/j ~A,B) and ground states |3^). As- 
sume that the atoms interact with the common vacuum and 
that transition dipole moments of atom A are parallel to the 
transition dipole moments of atom B. Due to this interaction, 
the process of spontaneous emission from two excited levels 
to the ground state take place in each individual atom but a 
direct transition between excited levels is not possible. More- 
over, the coupling between two atoms can be produced by the 
exchange of the photons. The evolution of atomic system can 
be described by the following master equation |0] 

^ = (L^+L« + L^> (ILl) 
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where for = A,B we have 

U'p =yi3 (2o^ipo^3 - o^3c/,'iP - P<30^i) + 

Y23 ( 2o^2P<^3 ~ ^i^nP ^ P<^3<^2 ) 

and 

L^«p = ri3 (2o^ipof3 - ofjo^ip - pofjo^i 



(11.2) 
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+ /£l23 [0f30f2+0f3a^2,P] 

In the equations (III.2l i and (III.3I I. O^^ is the transition operator 
from \kfj) to | ( /j = A,B ) and the coefficient 7^3 represents 
the single atom spontaneous - decay rate from the state | j) 
( j = 1,2 ) to the state |3). The coefficients Fjj, and Q.j^ are 
related to the coupling between two atoms and are the col- 
lective damping and the dipole - dipole interaction potential, 
respectively. As was shown in Ref. |7], in such atomic system 
there is also possible the radiative process in which atom A in 
the excited state \1a) loses its excitation which in turn excites 
atom B to the state |2g). This cross coupling between two 
atoms is sensitive to the orientation of the transition dipole 
moments of atoms, and in the present paper we study the 
model in which that coupling is absent. We also assume that 
the spontaneous - decay rates satisfy 

Yi2~Yi3=y (11.4) 
The remaining coefficients in equation (III.3I I can be written as 



(11.5) 



where j = 1 , 2 and R is the distance between atoms. A detailed 
form of the functions G/3(/?), Fjj,{R) depends on the geometry 
of the system [7], but in general, for R ^ °o 



andfor/?^0 



Gj3{R),Fj,{R)^0 



GMR)^l 



whereas the functions FjT,{R) diverge. 

The time evolution of the initial state p of the atomic system 
is given by the semi - group {7)},>o of completely - positive 
linear mappings acting on density matrices I113i1 . generated 
by +LP + L^^ . The properties of this semi - group cru- 
cially depend on the distance R between two atoms. As can 
be shown by a direct calculations, when the distance is large 
(compared to radiation wavelenght), the semi - group {7J}(>o 
is uniquely relaxing with the asymptotic state |3a) ® |3b). On 



the other hand, when R is small, Fn, r23 y and Hn, ^23 
are large, so we can use the approximation 
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In that case, the semi - group is not uniquely relaxing and 
asymptotic stationary states are nontrivial and depend on ini- 
tial conditions. 

We do not discuss details of the time evolution of the sys- 
tem, but we focus on the analysis of the asymptotic behaviour 
of the dynamics of atoms with small separation, when the 
conditions ( III.4b and ( III. 6b are satisfied. The master equation 
(III, lb can be used to obtain differential equations for the ma- 
trix elements of any state p. We consider matrix elements of 
p with respect to the basis of ® given by vectors 



1,2,3 



(11.7) 



taken in the lexicographic order. The equations for p/,„, l,m = 
I,..., 9 form a system of linear differential equations which 
can be solved by elementary methods. Using these solutions, 
after a long calculations, we obtain the explicit form of the 
asymptotic state p^s for any initial state p with the matrix ele- 
ments p/,„: 
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where 

x=\ (p22 + 2p33 + P44 + 2p77-2Rep24-4Rep37) 
o 

Z = ^ (P36 - P38 - P76 + P78) 

W = ^ (P26 + P28 + 2p39 - P46 " P48 ~ 2p79) (II-9) 
1 , 

y = g (P22 + p44 + 2p66 + 2p88-2Rep24-4Rep68) 
V = ^ (-P23 - P27 + P43 + P47 + 2p69 - 2p89) 

and 

t = \~2x~2y 

To get some insight into the process of creation of non - trivial 
asymptotic state p^s, it may be useful to consider the basis of 
collective states in C^, given by the doubly excited states 

kl)-|lA)®|lB), k2) = |2a) ® |2b) 

the ground state 

|^) = |3a)®|3b) 



3 



and generalized symmetric and antisymmetric Dicke states 
(see e.g. 01) 



1 



V2 



[\kA)(^\lB) + \lA)(^\kB 



(11.10) 



\aki) = -^[\kA) ^ \Ib) - \Ia) ® \kB)] 



where k,l — l,2,3;k < I. The states (III. 1 Oi l are entangled, but 
in contrast to the case of two - level atoms, they are not max- 
imally entangled. One can also check that the doubly excited 
states |ei), \e2) and the symmetric Dicke states \ski) decay to 
the ground state \g), whereas antisymmetric states jan) and 
|fl23) decouple from the environment and therefore are stable. 
Moreover, the state {an) is not stable, but it is asymptotically 
non - trivial. Notice that the collective states can be used to 
the direct characterization of the asymptotic behaviour of the 
system. In particular, the parameters x and y in dll.St are given 
by the populations in the antisymmetric states |fli3), |fl23) and 
\ai2)- 

x= ^ ((fli2|p|ai2) +2(ai3|p|ai3)) 

^((fll2|p|fll2)+2(fl23|p|fl23)) 

The remaining parameters can be computed in terms of the co- 
herences between the collective states. Since the populations 
(fli3|p|fli3) and (fl23|p|fl23) are stationary, the states which 
have the property of trapping the initial populations in {an) 
or |a23), create the non - trivial asymptotic state Pas with the 
stationary entanglement. On the other hand, the population 
in the state jo 12) is not stable, but it can be transformed into 
(ai3|p|fli3) and (a23|p|a23) in such a way that the values of the 
parameters x and y are fixed. The explicit examples of such 
behaviour of the system will be discussed in next section. 



in. GENERATION OF NPPT STATES 

To describe the process of the creation of correlation be- 
tween two atoms leading to their entanglement, we need the 
effective measure of mixed - states entanglement. As such a 
measure one usually takes entanglement of formation £'f(p) 
ifTsll . but in practice it is not known how to compute this mea- 
sure for pairs of li - level systems in the case when d > 2. 
A computable measure of entanglement proposed in |16] is 
based on the trace norm of the partial transposition p^^ of the 
state p. From the Peres - Horodecki criterion of separability 
|j8|,i9J, it follows that if p''^ is not positive, then p is not sepa- 
rable and one defines the negativity of the state p as 



A^(p) = 



,PT| 



- 1 



(III.l) 



A^(p) is equal to the absolute value of the sum of negative 
eigenvalues of p^^ and is an entanglement monotone, but it 
cannot detect bound entangled states 1,10.1 . Using the measure 



(IIII.ll l. one can check that generalized Dicke states are indeed 
not maximally entangled, since 

N{\s,i))=Ni\au))^^- 

In this section we study negativity of the asymptotic states 
(III.8I 1. For such initial states, where only the populations in 
antisymmetric Dicke states are non - zero, it is possible to ob- 
tain analytic expression for asymptotic negativity. By a direct 
calculation one shows that 



A^(pa 



y^4(x2+3;2)+f2_f 



(111.2) 



Notice that every non - trivial asymptotic state from this class 
is entangled. On the other hand, asymptotic negativity for 
initial states with non -zero coherences, can only be studied 
numerically. 



A. Pure separable initial states 

We start with pure states ( III. 7b . It is obvious that the initial 
states |1a) ® and \2a) \2b) decay to the ground state 
\g). On the other hand, the initial state |1a) CS) |3b) (atom A 
in the excited state and atom B in the ground state) has the 
population in the Dicke state jo 13) which is equal to j, thus 
for that state 



1 1 

x= - , t = — and y - 







and the asymptotic state is entangled with negativity 



A^(pas) = 



^/2~ 1 



(III.3) 



Similarly, the state \2a) ® \2>b) has the population 5 in the state 
|fl23) and also produces asymptotic state with the same value 
of entanglement. The same behaviour can be observed for the 
initial states |3a) (8" and \2>a)® \2b)- 

When two atoms are initially in different excited states i.e. 
we have the states \\a) ® \2b) or \2a) ®\^b), then the initial 
populations in the states jan) and |fl23) are equal to zero, but 
the population in the non - stable state \ai2) is non - zero and 
equals 5. During the evolution this population is transformed 
into the states |ai3) and \a2-i), in such a way that the asymp- 
totic state Pas satisfies 



(ai2|Pas|ai2) ^ 



and 



(fll3|Pas|fll3) 



(fl23|Pas|fl23) 



Thus the state Pas is also entangled, but its negativity is less 
then dlini l and equals (V6- 2)/8. 

Interesting examples of pure non - diagonal initial state are 
given by the following superpositions of states 1 1^) ® \2b) and 

^^ = cos(^|lA)®|2B)+sin(^|U)«)|3B), (^e [0,7t/2] 

(III.4) 
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FIG. 1 : Asymptotic negativity of ( IIII.4t as a function of (|). The values 
( -v/S — 2) /8 and ( \/2 — 1 ) /4 are indicated by dashed and dotted lines 
respectively. 



FIG. 2: Asymptotic negativity for the mixtures of states \ ® |2b) 
and \ Ia) ® \3b) (dotted line); |1a) ®\1b) and |1a) ® \3b) (dashed 
line); |1a)®|3b) and |2a)®|3b) (solid line). 



In that case, the initial populations 



1 



(flisl^Vjan) ^sin2(|), (fl23|^Jfl23) =0 



and 



(fll2|^'>pjfll2> 

are transformed into 



1 2. 
-COS 

2 ^ 



(fli3|pas|fli3) = :^sin^(|) + icos^(|) i (3-cos2(|)) 
z 4 o 

1 , 

(fl23|Pas|fl23) = ^ COS (|) 
(fll2|Pas(fll2| = 

In that case, the asymptotic states are parametrized by 

Jc= ^ (3-cos2(j)), 3;=icos-(|), f = ^, y = -isin2(j) 

and their negativity can be computed numerically. In FIG.l, 
we plot the asymptotic negativity as a function of (|). This fig- 
ure shows that for some values of the parameter (|), the super- 
position ( IIII.4I ) can have larger asymptotic negativity then that 
achieved by the initial state \1a) ^ |3b), which has the max- 
imal value of negativity produced by the dynamics for pure 
diagonal initial states. 



B. Some mixed separable initial states 

For the incoherent mixtures of pure states | Ia) (8) \2b) and 
I 1a) (8) |3b) i.e. the initial states 

p = p I U) ® |2b) (Ul ® (2b\ + (1 - /?) I U) «) |3b) (U I «) (3b| 

(III.5) 

the dynamics also produces entangled asymptotic states. One 
can check that their negativity is given by 



A^(pas) 



32 ^ ' 4 



(111.6) 



and observe that ( IIII.6I ) as the function of mixing parameter 
p is convex, so in contrast to the coherent superposition, the 
asymptotic negativity never exceeds ( 1111.3b . Analogous prop- 
erties of negativity can be found for another mixtures of two 
pure diagonal states (FIG. 2). 

To study the asymptotic entanglement produced for general 
class of mixed diagonal states, we consider two parameters 
to describe asymptotic states: their negativity and degree of 
mixture given by the linear entropy 



5L(p) = ^tr(p-p2) 



(in.7) 



As we show numerically, the set of all asymptotic states cor- 
responding to the diagonal initial states is represented on the 
entropy - negativity plane by the region Aas bounded by three 
curves (FIG. 3). The boundary curves can be found analyti- 
cally and are given by the following equations: the solid line 
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FIG. 3: The set Aas for the asymptotic states for all diagonal initial 
states. The solid line corresponds to the curve l llll.St ; the dotted line 
corresponds to l |III.9l l and the dashed line corresponds to dill. 10b . 
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on FIG. 3 is described by 



where 



Sl 



(III.8) 



64 



(8-.) 



for 5 e [3,8], the dotted hne is given by 



N- 



(111.9) 



fors G [1,2], and finally, the dashed line is described by equa- 
tions 



1 r- 1 
8 4 



(III. 10) 



Sl = 



64 



where [9/16,45/64]. 

Every asymptotic state corresponding to diagonal initial 
state is given by some point from the set Aas. In particular, on 
the curve ( IIII.8b lie asymptotic states produced from the mix- 
tures ofthe state |1a)<Xi|1b) and |1a)(8)|2b). Asymptotic states 
corresponding to the mixtures of |1a)(E)|1b) and |1a)® |3b) lie 
on the curve ( IIII.9I ). whereas on the curve jllLlOl i lie asymp- 
totic states obtained from the mixtures of states 1 1^) O |2b) 
and |1a)® |3b). 

An example of the state lying inside the set Aas is given by 
the asymptotic state generated from the initial state 



(III.ll) 



which is maximally mixed state of two qutrits. One can check 
that the asymptotic state is given by the formula (III. 81) with 
x = y— 1/12 and t = 2/3 . The corresponding values of linear 
entropy and negativity are equal to 9/16 and (3\/2 — 4)/12, 
respectively. As we see, even in that case, the incoherent 
process of spontaneous emission produces such correlations 
which entangle two atoms and diminish entropy of the sys- 
tem. 



C. Initial states with bound entanglement 



1 +fl 



a = 



and < fl < 1 . The states (IIII.12I ) have positive partial trans- 
position, but nevertheless are entangled jlQ}- Thek entangle- 
ment is bound and cannot be distilled ifllll . It can be checked 
by applying for example the realignment criterion of entan- 
glement 1 17, 18]. The criterion can be stated as follows: if the 
trace norm of realigned state R{p) is greater then 1, the the 
state p is entangled. We can also introduce the measure of en- 
tanglement based on this criterion. The so called realignment 
negativity Il9ll Nr{p) of the state p is defined by formula 

Nr{p) = \\R{p)\\-1 

This measure can detect some bound entangled states, but 
not all of them. In the case of states (IIII.12I I. the values of 
NR{pa) are contained in the interval (0, 0.0035) and the max- 
imal value is attained for a = | . 

Numerical analysis of the evolution of the initial states 
(IIII.12b indicates that their realignment negativity very rapidly 
goes to zero, but for the later times the states become entan- 
gled with positive values of negativity. Thus the dynamics 
studied in the paper has a remarkable property: bound entan- 
gled initial states ( |III.12t evolve into "free" entangled asymp- 
totic states. By a direct calculation one can show that the 
asymptotic states have the form (III.8b with parameters 
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The values of negativity of those asymptotic states are plotted 
in FIG. 4. 
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Consider now the following initial states 
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(III. 12) 



FIG. 4: Asymptotic negativity of the initial states l lIII.12b as a func- 
tion of parameter a. 



IV. CONCLUSIONS 

We have studied entanglement production in the system of 
two three - level atoms in the V configuration, coupled to the 
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common vacuum. In the case of small (compared to the ra- 
diation wavelength) separation between the atoms, the sys- 
tem has nontrivial asymptotic states which can be entangled 
even if the initial states were separable. Particular examples 
of such separable initial states are pure states in which the 
atoms are in different excited states. The process of the pho- 
ton exchange between the atoms, produces correlations that 
entangle two atoms. It is interesting that when we superpose 
such initial states, we can enlarge the amount the asymptotic 



entanglement. We have also characterized the entanglement 
of asymptotic states for mixed diagonal initial states. Using 
the description of mixed states in terms of degree of mixture 
and entanglement, we have found the region on the mixture - 
entanglement plane corresponding to such asymptotic states. 
We have also shown by considering a specific example that 
the dynamical evolution of that system brings bound entan- 
gled PPT states into "free" entangled NPPT asymptotic states. 
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